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Abstract 

We analyze the constraints imposed by Higgs vacuum stability on models with new fermions beyond 
the Standard Model. We focus on the phenomenology of Higgs couplings accessible at the Large Hadron 
Collider. New fermions that affect Higgs couplings lead to vacuum instability of the Higgs potential. Above 
the scale of vacuum instability, bosonic states must stabilize the potential, implying a cut-off to the pure 
fermion model. Conservatively tuning the models to produce the maximal cut-off for a given Higgs coupling 
effect, we show that observing a deviation in the Htt, f/-diphoton, or H-digluon coupling, larger than 20%, 
would require that new bosons exist in order to stabilize the Higgs potential below about 100 TeV. For generic 
parameter configurations, and unless the new fermions are made as light as they can possibly be given 
current experimental constraints, observing a 10% deviation in any of these couplings would suggest an 
instability cut-off below 10-100 TeY Similarly, if new bosons are absent up to a high scale, then a deviation 
in the Hbb or Htt coupling, larger than about 20%, should be accompanied by a sizable deviation in the 
Zbb orZTT couplings that can be conclusively tested with electroweak precision measurements at planned 
lepton colliders. 


1 Introduction 

The large hadron collider (LHC) Run-I gave us the Higgs hoson, hut the weak scale hierarchy problem does 
not seem closer to a solution than it did decades ago. This may change with new experimental information in 
Run-II, of which improved Higgs coupling measurements (T}|^ are a guaranteed outcome. A natural question 
to ask, is whether this Higgs particle that was found is the only one of its kind, namely, the only scalar particle 
up to very high energies. 

Indeed, proposed solutions to the hierarchy problem include new bosonic states beyond the Standard 
Model (SM). Examples are the scalar super-partners of the SM fields in supersymmetry and the bosonic 
resonances in composite Higgs models . In these examples, the scale at which the new bosonic states be¬ 
come dynamical marks the cut-off of the quadratic divergence in the quantum corrections to the Higgs mass. 

In this paper we show that measuring deviations in Higgs couplings at the LHC can establish the presence 
of new bosonic states, even if these bosons do not directly affect any Higgs coupling and are beyond reach of 
direct production. To show this, we proceed by elimination: we analyze the possibility that Higgs coupling 
modifications arise due to new fermionic states, without any new bosons up to a scale Ayy. Our task is then to 
derive an upper bound on Ayy. As we show, this upper bound is found from vacuum stability. 

To explain the logic, note that the only way to couple new fermions to the Higgs is through Yukawa cou¬ 
plings. New Yukawa interactions can certainly affect Higgs couplings to SM states, e.g. through new fermions 
running in the Hjy or HGG loop amplitudes, or mixing with the SM leptons or quarks at tree level. How¬ 
ever, as we shall show, the new Yukawa couplings must be sizable to generate a measurable deviation. Large 
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Yukawa couplings have a definite effect in the renormalization group evolution (RGB) of the Higgs-self quartic 
coupling, driving the quartic negative and leading to an instahility in the effective potential . To fix this 
instahility, at least within the domain of validity of a perturbative analysis, new hosonic states are neededj^ 
Vacuum stability has been invoked as a constraint on the SM effective theory in the past (for an early re¬ 
view see, e.g. [1^). A point that drew much attention in the days prior to the Higgs discovery and the precise 
measurement of the top quark mass, was the fact that measuring a heavy top quark or a light Higgs would have 
indirectly but robustly excluded the SM and required new physics to enter around a cut-off scale A that, for 
certain (then plausible) values of rut and m/j, could have been within reach of collider experiments such as 
the (then futuristic) LHC f^ -17 . Our paper can be thought of as an update circa-2015 of this logic. Today, 


having measured both the top and the Higgs masses to impressive accuracy (establishing that the SM Higgs 
potential is consistent with vacuum stability up to very high scales (Tsfl^ ), the missing crucial experimental 
information is the precise values of the Higgs couplings. 

By the end of the LHC 14 TeV program we expect uncertainties in the ballpark of 5-10% on Hbb, Htt, 
Hyy, HGG, Htt, HZZ, HWW with 300 fb“^ of data |^. Postponing more details to the body of the paper, 
our generic quantitative statement here is that a deviation at the ^(10%) level in any of the Hyy, HGG, or Htt 
couplings would imply that new bosonic states must stabilize the Higgs potential below a scale Auv ~ 100 TeV. 
An ^(10%) deviation in the Hbb or Htt coupling should be accompanied by a corresponding deviation in the 
Zbb or Ztt couplings at the permille level, well within the expected resolution of planned GigaZ machines 
such as the international linear collider (ILC) or circular electron-positron colliders in China or at CBRN; ruling 
out this corresponding Z-pole deviation would rule out pure fermion models, in addition, while we do not 
report a detailed analysis of this point here, applying our results to the HWW and HZZ couplings suggests 
strong vacuum stability constraints, strong enough to imply that observing a deviation in one of these channels 
at the LHC would rule out pure fermion models. 

Our results can be turned around to serve as generic prediction for the Higgs couplings in theories that do 
not contain new bosonic states up to very high scales, such as split supersymmetry |22||^ and its variants. 
According to our analysis, this general class of models predicts that Higgs coupling modifications will not be 
discoverable (or just very barely) at the LHC. This is not a trivial point because low-lying fermions, protected 
by chiral symmetries, could in principle be accommodated in these theories and couple to SM fields. 

Our results are relevant to the high-luminosity LHC as well as to a future lepton collider such as the ILC, 
the electron-positron mode of the future circular collider FCC-ee (formerly known as TLBP) and the circu¬ 
lar electron positron collider (CEPC), that promise percent and even sub-percent accuracy on Higgs coupling 
measurements 26 ^ . 

The body of this paper deals with the calculation of the cut-off scale Auv that could be inferred from mea¬ 
suring a deviation in the Higgs couplings to SM states. Since our derivation requires that we add only new 
fermions but no scalars or vector bosons, and since exotic chiral fermions are either ruled out already or, where 
they are not, require very large Yukawa couplings and so make our instability analysis trivial, we will only be 
dealing with new vector-like fermion representations. The fact that a sizable deviation in Higgs couplings in 
a pure-fermion theory renders the Higgs potential unstable was noted in several works (see also more 


general discussions in |36}j47) ). None of these works, however, have made the instability issue their primary 

and 


quantitative focus. The closest in spirit to our current analysis are Refs. |j^ and ||^ , where the instability cal¬ 
culation was done for the specific channels of Hyy and HGG. Here we expand on these works significantly by 
refining the instability calculation (using the one-loop, two-loop RGE-improved Coleman-Weinberg potential) 
and by generalizing to additional Higgs couplings. 


*An alternative logical possibility is that the Higgs scalar itself ceases to exist as a fundamental state above a cut-off scale Auv- This 
alternative is even more exciting but we are not aware of a framework in which this happens without new bosonic degrees of freedom 
(fundamental or composite) becoming dynamical close to the same scale. 


2 










The rest of this paper is organized as follows. In Sec. |^we consider Higgs couplings to SM fermions Hff. 
We begin hy analyzing the Hbb coupling in Sec. 2.1 and use this first example to illustrate our calculation of 


the Higgs effective potential and our definition of the vacuum instahility scale Auy. We then proceed with the 
Htt and Htt analyses in Secs. 2.2 and 2.3 We survey all of the new physics vector-like fermion representations 
that could mix with SM fermions at tree level, inducing Higgs-SM fermion coupling modifications. Analyzing 
together the Higgs and Z-hoson couplings, we show that the Hff and Zff coupling deviations are correlated. 
This leads to strong constraints on the Hbb and Htt couplings. In Sec. 2.2.1 we devote special attention to 
validating our analytical results with numerical calculations for the case of light vector-like fermions. In Sec.|^ 
we analyze the Hyy and HGG couplings. We conclude in Sec.|^ App.|^ provides details of our calculation of 
the Higgs effective potential. App.[^ contains details of the Higgs-fermion coupling analysis, surveying all of 
the relevant vector-like fermion representations that can mix with the SM states and presenting some useful 
formulae derived using effective field theory. In App.|^we give a rough, hut rather inclusive, assessment of 
the collider constraints on new physics fermions, based on the results of Run-I of the LHC. This survey of the 
current experimental constraints is important in particular for the Hff, HGG, and Hyy analyses. 


2 Higgs couplings to fermions 

Our task in this section is to consider the possible ways in which new vector-like fermions could mix with the 
SM fermions /, inducing 5rff^0 where 

[Shff - gfff) /gfff. ( 2 . 1 ) 

Only a few vector fermion representations can mix with the SM fermions. In App. we survey these repre¬ 
sentations, as well as define our notation for the SM and new physics fields. Higgs couplings and electroweak 
observables in some of these models have been studied extensively in the literature p^[^[50fj^ . LEP and 
LHC constraints imply that new charged fermions cannot be too light. In App. we summarize the existing 
collider constraints, finding that M > 600 GeV is a conservative lower limit on the mass scale of b and f quark 
partners in all cases of interest to us. This relatively heavy mass scale justifies the use of effective field theory 
(EFT) in analyzing the induced modifications to Higgs-quark couplings. For non-colored fermions the collider 
bounds are weaker, and we will devote some effort to extend the EFT analysis when dealing with Higgs-lepton 
couplings. 

Integration-out of heavy vector fermions up to non-renormalizable operators of dimension six is done in 
App. 1^ From this exercise one finds that the effective Higgs and Z-boson couplings to the SM fermions ex¬ 
hibit correlated modifications: large deviations in the effective Hff couplings imply sizable deviations in the 
Zff couplings. This leads to vacuum stability constraints that can be very relevant for the interpretation of 
upcoming Hbb and Htt data, as we show in detail below. In addition, for the top quark, the fact that the 
SM top Yukawa coupling is implies that large new physics effects are required to deform the Htt vertex 
appreciably, leading again to a strong vacuum stability constraint in the case of pure fermion models. 

2.1 Hbb 

To analyze the vacuum stability constraints on the effective Hbb coupling, and the interplay with precision 
Zbb data, we begin with a concrete model example. Consider the vector-like fermion representation 

Q(3,2)i, Q''(3,2)_i, D(3,l)_i, D''(3,l)i (2.2) 

6 6 3 3 
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with the potential 

= YQacH^Qd‘^ + YqD‘:ll^qD‘^ + YqdcH^QD^^ + YqcdH'^€Q‘^D + MqQ'^YM oDD^^ + cc. (2.3) 
This representation is denoted as rep’ DI in App.|^ 

Given a concrete example we can study the Higgs effective potential in the presence of the new fermions. 

In Fig. we plot i^ff as a function of the classical field he. We set Mq = Mo = 1 TeVi Yqcp, = 1.25, and 
= yqD'= = YQd<^ = 0, defined at the scale /r = 1 TeV. The details of the calculation of Kiff are given in App.|^ 
The vacuum instahility corresponds to the negative runaway of the effective potential, seen in Fig. (orange 
curve) to occur about an order of magnitude above the vector-like fermion mass scale. 

\/eff/(10 TeV)^ 
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Figure 

of Eqs. 


1 : The one-loop Higgs effective potential, shown for the SM (blue) and in the presence of the vector-like fermions 
2.2[|2.3 (orange). We set Mq — Md — 1 TeV, Yqcp — 1.25, and igu' — Yqoc — Yqd^ — 0, defined at the scale p = 1 TeV 


The vacuum instability depicted in Fig. can be understood as due to the negative RGB running of the 
Higgs quartic coupling in the presence of extra fermions. In Fig. we plot the Higgs quartic coupling A for the 
same model example of Fig.[^ In the plot, the blue line shows the SM running of A and the orange line gives A 
with the new fermions included. The jump at /i = 1 TeV is due to the threshold correction in matching the SM 
EFT below the vector mass scale to the full theory above it. 

Figs. [T]j^ demonstrate that the field content of Eq. j2.2) cannot be considered as a complete theory, if im¬ 
plemented with a large Yukawa coupling at the vector-like fermion scale. Other fields, coupled to the Higgs, 
must be added in order to eliminate the runaway of the effective potential. Adding more fermions would make 
the instability worse, so the new fields need to be bosons. 

A precise determination of the mass scale at which the new bosons need to be introduced, given an ap¬ 
parent runaway in the effective potential as in Eig.[^ is beyond the scope of this paper. Instead, we define a 
rough criterion for the scale Ayy at which the effective potential should be stabilized, through the condition 
Aeff(Auv) = —0.07 where Aeg parametrizes the effective potential at large he values, where we can set all di¬ 
mensional parameters to zero besides from the Higgs field itself: I^ff —> ^ value Aeff = —0.07 indicates 

roughly the onset of vacuum instability; for more negative Aeff, the time scale for tunneling from he = 246 GeV 
to a remote vacuum with he > Auv becomes exponentially shorter than the age of the Universe [53]pj In Fig.|^ 
we plot Aeff for the same example of Figs. EB 

rough estimate of the tunneling probability through true vacuum bubbles of nucleation size f/A is p ~ {A/Hoy where 
S(A) ^rid Ho ~ GeV is the ffubble constant. Setting p = l gives A(A) —0.065 (l -I- 0.02logj„ (A/1 TeV)) \ 
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Figure 2: The Higgs quartic coupling in the SM (hlue) and with the model example of Fig. (orange). 



Figure 3: The effective Higgs quartic coupling in the SM (blue) and with the model example of Fig. (orange). The 
horizontal green line marks Aetf = —0.07 indicating a rough criterion for vacuum instability. 
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The scale Auv> or any similar characteristic scale parametrizing the potential runaway, depends on the 
size of the Yukawa couplings. For Yqco slightly larger than unity, as in the case depicted in Figs. [T]j^ we find 
Auv about an order of magnitude above the vector-like fermion mass scale. Larger values of Yqcd within the 
perturbative window (| Y| < 47 i:/vC^, with Jf = 4Aic = 12 in the current case) change Auv by an order one 
amount. Similarly, turning on, in addition to some of the other Yukawa couplings in the model, also has a 

modest order one effect. To substantiate this point, in Fig.|^we plot the ratio Auv/Afi, where Mi is the lightest 
vector-like fermion mass eigenstate, vs. the Yukawa coupling Fqc^,. In obtaining the green (blue) smooth fines, 
we set Y = Yqc^, all other Yukawa couplings to zero, and Mq = Mjj = 1 TeV (5 TeV). In the dotted lines we 
repeat the same values for Mq,Md, but turn on the additional Yukawa coupling Yq^c = 1. We have verified that 
varying Mq 7 ^ Mu does not change the results. 



Figure 4: Auv/Mi vs. the Yukawa coupling Y —Yqc u. The model is the vector-like quarks of Eqs. j2.2|2.3 . Smooth (dotted) 
lines show the results fixing Yq^c, an additional Yukawa coupling in the model, to Yqqc — 0 {Yq^c = 1). Blue (green) lines 
are for Mq = Mn = 1 TeV (5 TeV). 


The bottom line is that allowing any of the Yukawa couplings in Eq. 2.3 to become larger than ~ 1.5 would 


render the vector-like fermion model inconsistent right above the new fermion mass scale. However, it is also 
apparent in Fig. that if the Yukawa coupling is made sufficiently small, then the instability can be pushed 
up to high energies. This decoupling of Auv, that occurs as the Yukawa couplings are made small, is important 
because it means that for certain parameter configurations the vector-like fermion model can form a consistent 
EFT up to a high scale, in which case the vacuum stability argument looses phenomenological relevance. 

The next step is to exhibit the relation with Hbb and Zbb measurements. Integrating out the heavy fields 


we find (see 1 32 ^ for related analysis) 


5rb 




,icj> 


c/Zrub 


(2.4) 


where 0 is a complex phase (to be explained shortly) and 5gvb and 5gAb are the modification to the vector 
and axial Zbb couplings, respectively, defined in the usual way (see, e.g. |j^) and constraint by LEP data to 
be in the ballpark of a few percent . In App. (see Tab. we give a derivation of Eq. l |2.4| by integrating 
out the heavy states in an expansion to leading order in Y^v^/M^, where Y is any of the Yukawa couplings 
in the problem and M = Mq^d- At leading order in Y^v^/M^, as done in App. the phase 0 is given by 
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(j) = arg^iQC£)i^£,ci^rfcy*MQM^j. However, and usefully for our purpose, Eq. 2.4' is actually valid to good 
accuracy even when the vector-like states are light and allowing the Yukawa couplings in the vector-like sector 


to be large, namely YqjjcV^/M^, Yqcj^v^/M^ = To see this, note that one could also derive Eq. i2.4' ex¬ 


panding in the mixing couplings Yqjjc, Yq^c , and in the mostly-SM bottom Yukawa coupling yh, but keeping all 
orders in Yqco and The couplings j/fo, igds are guaranteed to be good expansion parameters because 
of the smallness of the bottom quark mass (combined with the fact that 5 is constrained to be smaller than 
unity by existing LHC data) and by the experimental constraints on 5 gAb.vb ■ Following this route and allowing 
Yqjjc v^/M^, YqcjjV^/M^ = &{!), we hnd that Eq. |2.4| remains valid as is, and the only modification that needs 


to be done is a generalization of the phase cj) = arg ^ 


YqC pY^pC Yq^C MqMd^ 


where 


is the determinant of the 


3x3 mass matrix given by Eq. j2.3| and \M 2 ^ 2 \ = MqMo — YqcjjYqoc v'^12 is the determinant of the vector-like 
2x2 sub-matrix. In Sec. 2.2.1[ devoted to vector-like leptons but dealing essentially with the same formula, we 
demonstrate the validity of Eq. |2.4| for very light vector-like states with a numerical analysis. 


Eq. 2.4 shows that a large deviation in Hbb requires some corresponding deviation in the Zbb couplings. 
This correlation can be relaxed, but only at the cost of a sizable Yukawa coupling, Yqcd in our example. To make 
a quantitative estimate, note that for |5rb| > 0.2 or so, we can neglect the —25gAd contribution on the RHS of 
Eq. \2.4\ . Setting 5gAd\ ^ 10“^, within the expected resolution of future lepton collider experiments, we 

have Idrfal 0.1|iQC£)| j. Thus, given our analysis above summarized in Fig.j^ vacuum 

stability can make a powerful discriminator for a consistent explanation of Higgs couplings deviations, once 
the Zhh couplings have been determined to about a permille accuracy. 

In Fig.j^we demonstrate our results quantitatively. On the x- and y-axes we plot 5gvb and 5gAb, respec¬ 
tively. Inside the orange-shaded region, the coupling \ Yqcj)\ needs to be larger than 1 in order to induce an 
Hbb deviation of brj, > 0.2. Inside the gray-shaded region, \Yqco\ needs to be larger than 1.5 to achieve the 
same br^. To compare with current Zhh data, the blue contour shows the 95%CL allowed range for 5gvb,Ab’ 
using the reduced covariance matrix from Ref. that includes the anomalous LEP result for The green 
contour shows the 95%CL allowed region, obtained, ignoring correlations, using the measurements for Rb and 
Ab as taken from |^, and omitting the measurement ofA^^. 

We learn that the current Zhh coupling measurements still allow a signihcant Hbb deviation consistent 
with vacuum stability. Indeed, in vector-like quark interpretations of the A^^ anomaly |^, an Hbb deviation 
would be generically expected]^ Improved experimental determination of the Zbb couplings, in the ballpark 
of planned future experiments, can sharpen these results and test conclusively a fermionic interpretation. To 
compare with future Zbb data, the dashed black contour shows an estimate of the 95%CL allowed region fol¬ 
lowing future measurements at the ILC or another experiment of comparable precision, for which we adopt 
the SM central values for Rb and Ab and assume o-r^^ = 0.00014, (Taj, = 0.001 (to be compared with the 


current aRf^ ^ 0.0007 and a a,, ^ 0.02 from LEP data 155]). The right panel shows an expanded version of the 


left. The fact that the dashed black error ellipse is completely contained in the orange shaded region, signaling 
low scale vacuum instability, shows the potential power of the vacuum stability argument to constrain Higgs 
coupling deviations in conjunction with improved Zhh data. 

Before we conclude this section, we pause to generalize the results to other vector-like fermion represen¬ 
tations that can mix with the SM b quark. These representations are collected in App.|^ Tab.[^ For all of these 
cases, the vacuum stability analysis proceeds in a similar manner, leading to results similar to what we have 
shown in Fig.j^ Considering the Zbb couplings, models DII — DIV in Tab.[^possess the same basic structure 
as in Eqs. f2.2p.3| and produce relations between the Hbb and Zbb effective couplings, that are analogous to 
Eq. (2.4 1 and that were illustrated in Fig.[^ For completeness, we summarize the Hbb and Zbb relations for all 


^Note, however, that for rrih ^ 125 GeV vacuum stahility analysis indicates that the models of exhibit instability on scales 

much lower than the scale of gauge coupling unihcation. 
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Figure 5: Left panel: blue solid, green solid and black dasbed contours mark the 95%CL allowed region in the 5gvb,5gAt 
plane, considering current LEP Zbb data including LEP data omitting and future precision measurements from 

an experiment such as the ILC, respectively. Inside the orange-shaded (gray-shaded) region, the coupling needs to 
be larger than 1 (1.5) in order to induce an Hbb deviation of 5 ri, > 0.2. Right panel: close-up of the same results. 


of these models in Tab.|^ In contrast, model DV provides only one vector-like fermion state that mixes with the 
SM right-handed b quark, with no corresponding partner for the left-handed SM quark. As a result, the Hbb 
coupling and the Zbb coupling are directly correlated in model DV, predicting 5r},^ —25gvb- Using the cur¬ 
rent LEP constraints, model DV is limited to induce at most |5rfo| < 0.1 at 95%CL. This discussion makes clear 
why we do not devote attention, in the context of the Hbb coupling, to vector-like fermion representations 
in which only either the left-handed SM quark or the right-handed SM quark mixes with new fermions, hut 
not both. Examples for such models are given by simply deleting, e.g., either the pair Q,Q‘^ or the pair D,D^ 
in Eq. f2.2) , and similarly for the other cases DII — DIV. For all such vector-like fermion models, the Hbb 
deviation is directly tied to a corresponding Zbb deviation, implying that Idr^l > 0.1 is ruled out by existing 
precision Z-pole data. 

Finally we comment on loop corrections to Eq. \2A\ , that was derived at tree level. The dominant loop effect 
is due to the renormalization of the bottom quark Yukawa coupling from the scale to the scale mu, relevant 
for h ^ bb decay. This effect is captured by using the running MS bottom quark mass mb{mh) ^ 3 GeV in 
Eq. \2A\ . Another effect is the RGE evolution, within the SM EFT, of the nonrenormalizable operators produced 
by integrating out the vector-like quarks at their mass threshold down to the scale ruh [^60 . This running 
corrects Eq. \2A\ at the 10% level, in the direction of suppressing the right-hand side of Eq. \2A) compared to 
its value when using | iQ^ul as given at the vector-like quark mass threshold. 

To summarize: if a pure-fermion model produces a deviation > 20% in Hbb, it should also produce an 
observable deviation in Zbb. If the former if found without the latter, then bosonic states should exist in the 
spectrum not far above the fermion mass scale. 
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2.2 Hzz 


Our analysis of Hzz follows closely the Hbb discussion of the previous section. Again, we work out the details 
for one representative model example, assuming the field content 

L[l,2)_i, L"(l,2)i, E{1, l)_i, £"(1, l)i (2.5) 


and the potential 


= YiecH^Le'^ + Yi + Yle^H^LE^^ + YlceH^ eHE + MlV eH + MeEE‘= + cc. 


E 1 TJC . 


E 


T^tC t 


T^tC 


( 2 . 6 ) 


We consider the Higgs vacuum stability first. In Fig. we plot the ratio Auv/Mi, where Mi is the lightest 
vector-like fermion mass eigenstate, vs. the Yukawa coupling Yl^e- In obtaining the green (blue) smooth lines, 
we set Y = Yl^e, all other Yukawa couplings zero, and Ml = Me = 1 TeV (5 TeV). In the dotted lines we repeat 
the same values for Ml, Me, but turn on = 1. Varying Ml ^ Me does not change the results. We see 
that allowing any of the Yukawa couplings in Eq. |2.6| to become larger than ~ 2 would render the vector-like 


fermion model inconsistent right above the new fermion mass scale. 



Figure 6: Ayv/Mi vs. the Yukawa coupling Y = Yicp. The model is the vector-like leptons of Eqs. 2.5 


2.6 . Smooth 


(dotted) lines show the results fixing Ylec , an additional Yukawa coupling in the model, to = 0 (= 1). Blue (green) 
lines are for Ml — Me— I TeV (5 TeV). 


We next explore the relation with precision Zzz measurements. Integrating out the heavy fields, we hnd 
(see App.|^ Tab.[^ 


SVt: 


+ 


2\Ylce\v 

'jT.m.x 




At 


,i(p 


(2.7) 


The Zzz coupling deviations 5gw and SgAr are constrained by LEP data to be in the ballpark of a few per- 


mille . In analogy with our results in Sec. 2.1 Eq. 1 2.7 is valid to good accuracy even when the vector-like 
states are light and the Yukawa couplings Ylce and Yj,£c are the only subtlety being the proper definition 
of the phase (p. When the vector-like states are heavy we have (p = arg(^Yj,c£Y/£cl^ecy*M2M^j, while if one 
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allows light vector-like states then (p needs to he generalized in exact analogy with the discussion in Sec. 2.1 
In the rest of this section we first use the analytical result in Eq. f2.7) , and subsequently turn to a numerical 
verification for very light vector-like states. 

In Fig.j^we illustrate the vacuum stability constraint vs. Ztt data, computed using Eq. l |2.7| , inthe5gyT,5g,4 
plane. Inside the orange-shaded (gray-shaded) region, the coupling | needs to be larger than 1.5 (2) in or¬ 
der to induce an Htt deviation of > 0.2. The green solid contour shows the 95%CL allowed range for 
5gvT:,Ar, using the measurements for and at LEP To compare the result with future Ztt data, 
the black dashed contour shows an estimate of the 95%CL allowed region from future measurements at the 
ILC/GigaZ or another experiment of comparable precision, for which we adopt the SM central values for 
and At; and assume = 0.004, aA^ = 0.001 (to be compared with the current 0.045 and 0.004 

from LEP data ). The estimate of future ctr^ at the ILC/GigaZ is taken from Ref. . There is no discussion 

of (JAt at future lepton colliders in the literature; here we simply assume a value of 1/3 of the current system¬ 
atic error of A^. We expect that future lepton collider measurements would reduce the statistical error to be 
negligible compared to the systematic error. Ghanging the phase of 5 r-r does not appreciably affect the results. 



Figure 7: Green and black contours mark the 95%CL allowed region in the 5 gi/r, dg^x plane, considering current LEPZtt 
data, and future precision measurements from an experiment such as the ILC, respectively. Inside the orange-shaded 
(gray-shaded) region, the coupling | I/cbI needs to be larger than 1.5 (2) in order to induce an Hzz deviation of > 0.2, 
indicating vacuum instability. 


To summarize: a new physics Yukawa coupling | Y| > 1.1 implies that the Higgs effective potential goes 
unstable on a scale about three orders of magnitude (or less) above the vector-like lepton masses; | Y| > 2 leads 
to instability about one order of magnitude above the vector-like lepton scale. Absence of signals in future 
Z-pole precision measurements would exclude the pure fermion model from generating Idr^l > 0.2. 
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2.2.1 Numerical verification for very light vector-like states 

Vector-like leptons, especially if they only mix with the t and not with the electron or muon, are experimentally 
allowed to he quite light (see App.[^. We therefore turn to examine numerically the validity of Eq. \2.1\ for very 
light vector-like states, with the goal of validating Fig.[^ To do this, we numerically scan over the parameter 
space of the model in Eq. j2.6| . To make the scan efficient, we reparameterise the model by writing the 3x3 
lepton mass matrix M as 




-if = [rL-E-)ji 


E<^+ I , ^ = VDU, 

y — gi0„A7g!Y„A3g!^„A2g!<>A3^ JJ — giPuMgiruMgiSuM ^ 


( 2 . 8 ) 


where D = diag(mT,Mi,M 2 ) is a positive-definite eigenvalue matrix (with Mz > Mi) and V and U are unitary 
diagonalization matrices that we express using a convenient basis of the Gell-Mann matrices A, with real pa¬ 
rameters Pv^u,dv,u,Tv.u>(t>- This parametrization is convenient and minimal, making use of the [SU(2)xU(l)]2 
global symmetry of the gauge-kinetic terms in the action. It allows us to enter physical mass eigenvalues as 
input, and because the physical Ztt couplings are given simply by dgz, = | |i§il^ and 5gR = - | It/ial^, allows 
us to skip phenomenologically unacceptable model points from the outset. Moving back and forth between 


the basis we use for the numerical scan, Eq. 1 2.8 , and the basis of Eq. j2.6| is straightforward, and we do it to 
evaluate the size of the Yukawa couplings as defined in Eq. j2.6| in order to make the connection with vacuum 
stability. In the scan we omit CP-violating phases, setting (j) = Yv = Tu = 0- Besides from this simplification, we 
scan over the physical heavy mass eigenvalues Mi and Mz and over the angles 0^,^ and Eor each point 
in the scan, we evaluate and Sgw^r numerically, and, in addition, find the value of the maximal Yukawa 
coupling as given in the basis of Eq. j2.6| . 

The results of the scan are given in Fig.[^ where we repeat the calculation of Fig.j^and superimpose as blue 
points all of the scan points in which > 0.2, the maximal Yukawa coupling is smaller than 1.5, and the Ztt 
couplings are consistent with current data at the 95%CL. Thus, the scan points in Fig. [^should complement 
the orange shaded region, within the 95%CL contour. As can be seen from the plot, the results of the analytical 
study are confirmed by the numerical scan, and apply also for very light vector-like lepton states, here allowed 
to be as light as 200 GeVi with sizable Yukawa couplings. 

As final comments, we note that the occurrence of the scan points in the top region in Fig. is due to the 
exact relations 5gv = ^ and 5gA = ^ (if^il^ + llTisI^)- Other vector-like lepton representations 

that mix with the tau satisfy analogous expressions to Eq. \2.7\ , that can be found in App.|^ Tab.[^ Requesting a 
large 5 and imposing an upper bound to the Yukawa couplings, all of these representations are also restricted 
to specific regions in the 5 gw, 5 gAT plane: rep’ LI and LI 11 populate the top region shown by the blue points 
in Fig.[^ rep’ LI I and LIV populate the middle-right region, with > 0 and dgyiT around zero; and rep’ LV 

populates the middle-left with 0 > 5gw and SgAr around zero. In all cases, the shaded orange or gray regions 
of Fig. [^indicate large Yukawa couplings and, therefore, vacuum instability. 

The numerical analysis that we report here applies with obvious changes to the Hbb case as well, and 


demonstrates the validity of Eq. \2A\ on which we based the discussion in Sec. 2.1 


2.3 Htt 

The top quark Yukawa coupling in the SM is Gombined with the existing collider limits on fermion top 
partners, this means that large Yukawa couplings are necessary to give appreciable corrections to the Htt 
vertex through mixing with vector-like fermions, implying a strong vacuum stability constraint. 
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«5gvr 

Figure 8: Same as in Fig.[^ but showing in blue points tbe results of the scan described in the text to validate the EFT anal¬ 
ysis for low-mass vector-like states, where we impose |7| < 1.5. We scan on the parameter space of the model, allowing 
the lightest vector-like fermion state to be as light as 200 GeV. 


Consider the vector-like fermion representation 

Q(3,2)i, Q"(3,2)_i, U{3,1)2, U^{3,1)_2 

6 6 3 3 


(2.9) 


with the potential 

yNP = YQucH'^eQu‘^ + Yqu^H'^€qU‘' -E YqucH'^eQU'' -E YqcuH'^Q'^U+ MqQ^ eQ‘^ + MuUU‘^ -E cc. (2.10) 


This is rep’ UI inApp.j^ Due to LHC constraints, the vector masses Mq, Mu must he larger than about 700 GeV. 
Integrating out the heavy fields, we find 


2 f \YquA^ \YquA^ YqcuYqucY^u'^ A 

1^2|My|2 + 2|MQ|2^ MqMu ]’ 


( 2 . 11 ) 


which implies 


\5rt\ < 0.06 


Yquc (700 GeV) 


Yquc {700 GeV) 
Mq 


2 

+ 2\Yqcu\ 


Yqu'^ Yqw^ (700 GeV)2 
MqMu 


( 2 . 12 ) 


We can make a conservative estimate of the maximum deviation in Htt, consistent with a pure fermion 
model, hy letting all of the new Yukawa couplings he ~ 1, tuning their phases to interfere constructively as in 
Eq. I 2.12| , and at the same time allowing Mq,Mu ~ 700 GeV This gives \5rt \ < 0.25 and, judging from Fig.|^ 
that represents similar RGE, implies Auv ^100 TeV 

Note that the representation above gives a non-minimal example in terms of the field content. In contrast 
to the Hbb and Htt cases, where this possibility is precluded by the stringent Z-pole constraints, here we 
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could omit one of the pairs of vector-like fermions (either the Q, Q‘^ or the U, U‘^ pair) and still produce a po¬ 
tentially sizable deviation in the Htt vertejQ The result would correspond to Eq. 2.12 with one of Yqu^ or 
Yquc set to zero. This does not change the conclusion. Finally, all of the vector-like fermion models listed in 
App.j^yield similar results to the example we analyzed here. 


3 Higgs couplings to massless gauge bosons 


We now study the implications of vacuum stability constraints on HGG and Hyy coupling modifications in¬ 
duced by vector-like fermions. In this section, we ignore the mixing of the new fermions with the SM fermions, 
which was analyzed in Sec. Vacuum stability constraints in vector-like fermion models modifying Nyy and 
HGG were considered in and ||^. Our results here refine those analyses using the two-loop RG improved 
Coleman-Weinberg potential, as described in App.^ 

The leading-log contribution of the new fermions to the HGG and Hyy amplitudes can be derived from 
the Higgs low energy theorem 1 63p4 . For a fermion carrying electric charge Q, and transforming under SU(3) c 
in a representation of dimension D with Dynkin index T, we have 


5rr 


5 re 


4Q2D fdlog\\M\\'\ 


51ogi^ J’ 
piogllMlh 
I d\ogv )’ 


(3.1) 

(3.2) 


where ||M|| is the absolute value of the determinant of the fermion mass matrix M and —6.5 is the SM 

Hyy amplitude. Our convention is such that 7(3) = 1/2. 

The minimal building block of interest for our vector-like fermion representation is im 

i(D,1)q, (3.3) 


with the mass matrix for charge Q states 






(3.4) 


corresponding to the potential given in App.|^ Eq. i A.2 1 . In what follows we ignore the physical complex phase 


arg(m^m* and assume a basis where all of the parameters are real. We further choose m^,m^ to be 
positive. Denoting the mass eigenvalues as Mi and M 2 with 0 < Mi < M 2 , the log-determinant derivative is 
given by 


^log||M|| ^ 

d\ogv J 


YY‘^1’^ 

M1M2 


sign(|M|). 


(3.5) 


2 

Note that |M| = — y FF"^, so that a negative sign for |M| in our basis requires large Yukawa and small 

vector-like fermion masses. For >174 GeV, the product FF'^ must be larger than unity to have negative 

|M|, so that such configurations are automatically associated with low-scale vacuum instability. 

''Ref. |62| derived indirect constraints on the Zf f couplings that read ISgy,|, \5gAt\ S 0.1. Thus a partial vector-like fermion model 
with either Q,Q'' or U, U’^, hut not hoth, could in principle induce |5r,| ~ 0.2, compatible with electroweak precision data, hut at the 
cost of vacuum instahility as shown in the text. 

^For real color representations and Q = 0 or Q = we could set or !/> = t/)'. 

®We could promote x and /“'to SU (2) w triplets, which would add more fields to the RGB and would not change our conclusions. 
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Before proceeding to the numerical results, we comment on the effects of NLO corrections. The first NLO 
effect involves two-loop diagrams that contribute to the Wilson coefficients of the effective dimension six 
\H\^GG and \H\^FF operators at the vector-like fermion mass threshold. These corrections maybe sizable, of 
order tens of percent, especially for the HGG case (see, e.g. |[^ for the analogous effect in the integration out 
of the SM top quark). However, note that as long as we restrict our discussion to fermions in the fundamental 
representation of color, the dominant effect is a multiplicative factor that acts on the SM top quark amplitude 
and new physics LO contribution alike, and so drops out in the relative correction to the HGG vertex that we 
discuss here[^ Once we consider higher color representations (in the second part of Sec. |3.l| below), the NLO 
K-factor ceases to be a common multiplicative effect for the SM and new physics contributions, and would 
alter our results to some extent. The second NLO effect pertains to the RGB running of the nonrenormalizable 
operators from the vector-like fermion mass scale down to m/j in the SM EFT |^. We have checked that the 
corresponding correction is limited to a few percent, and we omit it in what follows. 


3.1 HGG 

In Fig. 1^ we show the maximal deviation Src as a function of the vacuum instability scale Ayy, assuming a 
vector-like pair of fermions in the fundamental representation of SU(3)c. The smooth purple, red, and orange 
curves correspond to setting the mass of the lighter colored fermion Mi to 0.5, 0.7 and 1 TeV; respectively. 

In producing Fig.|^we use the following method. Typically, the most conservative vacuum stability con¬ 
straint (largest Auv) is achieved when the new Yukawa couplings are as small as possible for fixed 5 re and Mi. 
Considering dre > 0, we see that Eqs. l |3.2| and |3.5| require YF^ < 0. The smallest consistent choice for |F| 
and I F^l, that we employ in Fig.[^ then corresponds to F —F^, with Mi M 2 and as small as is allowed by 
direct collider searches for colored fermions. Considering Stq <0, we need FF*^ > 0. The most conservative 
configuration for the vacuum stability analysis is different than that in the 5rc>0 case, due to the inequality 


(M2-Mi)2> 


{Y+Y‘^fv^ 


(3.6) 


which can be derived from Eq. f3.4| . Due to Eq. j3.6| we cannot tune F F*^ and Mi M 2 at the same time. 
Instead, the most conservativ e co nfiguration for the stability analysis in this case, that we employ in Fig. 


This explains the asymmetry o 




corresponds to saturating Eq. 13.6 1 with F F'^, | F| M 2 = Mi 


5 re I / gre ro~ 
2T \ 2T V 


Srr, t 
2T J 


'Luv for fixed in Fig.j^between positive and negative values of drc. 

We comment that the intuition by which larger Auv corresponds to smaller Yukawa couplings, that we used 
to fix the model parameters in Fig.|^ holds well as long as the instability threshold occurs sufficiently far from 
the vector-like fermion mass threshold. When the instability occurs immediately above the vector-like fermion 
mass scale, we find in some cases that larger Yukawa couplings can lead to a slightly higher instability scale, 
due to threshold effects. We stress, however, that these effects only become relevant when the instability scale 
is very low in the first place, Auv < 10 TeV. We comment about these effects further in the next section. 

In producing Fig.[^ we restrict our calculation to the region of parameter space where perturbation theory 
is under reasonable control by imposing |F|, |F^| < -^=. For fermions with D = 3 we thus impose |F|,|F^| <3.6. 
For most of the curves shown in Fig. all of the couplings remain perturbative at all RGB scales up to the 
vacuum instability scale Auv- An exception occurs for the M\ = 1 TeV example with large negative Stg < —0.1. 
Here we find that F and F^ run large with increasing RGB scale, and cross the perturbativity threshold defined 
above before our nominal criterion for vacuum instability applies. When this happens, we redefine Auv as 
the scale where the perturbativity threshold was crossed, which leads to the kink in the plot. Of course, this 


^For an NNLO computation in a related model (partial vector-like quark representations), confirming these statements, see 
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Figure 9: The maximum HGG coupling deviation vs. the corresponding instahility scale Auv, obtained hy adding vector¬ 
like fermions in the 3 representation of SU(3)c. Purple, red, and orange lines correspond to fixing Mi = 0.5,0.7 and 1 Tey 
respectively. The dashed line, that is only visible for the Mi —I TeV case at large negative 5 rg , demonstrates the sensitivity 
of the result in this parameter region to our criterion for perturbativity; see text for details. 


procedure is somewhat arbitrary; the main lesson is that for such large negative 5re, the model runs strongly- 
coupled quickly, making our perturbative estimates less reliable. To highlight the sensitivity of our results in 
this parameter region to changing the perturbative prescription, we superimpose as a dashed line the result 
obtained when modifying the numerical perturbativity threshold by 20%. 

Higher color representations have a larger Dynkin index and thus reduce the size of the Yukawa couplings 
needed for a given 5 tq, relaxing the vacuum stability constraint. At the same time, collider constraints on Mi 
become stronger (see App.[^, balancing the effect to some extent. 

In Figure [To] we show the vacuum stability bound for vector-like fermions in the D = 6 [T = 5/2, left panel) 
and D = 8 [T = 3, right panel) representations, using Eq. 13.3 1 with Q = 1, and taking Mi = 1.2,1.5 and 2 TeV. 
We do not consider representations of D > 8, as these cause the SU(3)c gauge coupling gs to run strong quickly 
and hit a Landau pole on scales very close to the vector-like fermion mass. 

Once again, in Fig. 10 we restrict our calculation to the region of parameter space where | Y|, | F^| < For 
d re > 0, all couplings are perturbative according to this criterion up to the vacuum instability scale. For large 
negative 5 rc , however, in both the D = 6 and D = 8 examples, the Yukawa couplings run strong and cross our 
perturbative reliability criterion at a scale lower than the vacuum instability. We indicate where this happens 
by showing, in dashed lines, the results obtained while modifying the perturbativity criterion by 20%. 


3.2 Hyr 

The largest effect in Nyy at a given stability cut-off Auv is obtained with color singlet fermions. This is because 
LHC searches put strong mass constraints on exotic colored fermions, whereas color singlets are still allowed 
to b e qu ite light (see App.[^. For example, an exotic color octet vector-like representation would give D = 8 in 
Eq. 3.l| , but this would come at the price of a large mass suppression. Mi > 1 TeV implying | | ~ 0-06 

even with large Yukawa couplings Y, ~ 1. In contrast, vector-like leptons with D = 1 are still allowed with 


Ml ~ 200 GeV and 


^glogllMlAl 


1 for the same F, ~ 1. The increase in 


pIogllMlA 

V S\ogv j 


d\ogu j 

for the decrease in D. Therefore in what follows we focus on color singlet representations. 


more than compensates 
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Figure 10: The maximum if GG coupling deviation vs. the corresponding instability scale Ay v, obtained by adding vector¬ 
like fermions in the 6 representation of SU(3) c (left) and in the 8 representation of SU{3)c (right). Purple, red, and orange 
lines correspond to fixing Mi = 1.2,1.5 and 2 TeV, respectively. The dashed lines correspond to reducing the perturbativity 
constraint by 20%. 


Given a value for Q and having set D = 1, we can calculate the maximal Ayy for a given 5 ry in analogy with 
the 5 re case. Typically, the maximal instahility scale for 5 < 0 is obtained for F — Y‘^ and Mi M 2 , while for 


5 ry> Owe have it for F F'^ and 


M 2 — Ml ^1 + X + (2 -P x)^ 


withx = 


5ry/4Q^D\. We comment again 


SM^'r 


that these relations are home out in our calculation as long as the instahility scale is sufficiently far from the 
vector-like fermion mass threshold, and imply simply that larger Yukawa couplings trigger an earlier vacuum 
instahility. In the extreme case in which the instability scale occurs very close to the vector-like mass, we find 
that threshold effects can conspire to make a slightly higher instability scale pair with slightly larger Yukawa 
couplings. 

In the left panel of Fig. 
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we plot 5 ry vs. Auv for a vector-like lepton representation, given by Eq. l|3.3| with 


D = 1 and Q = 1. Purple, red, and orange lines correspond to fixing Mi = 200,400 and 600 GeV; respectively. In 


the right panel we set Q = 2 and plot the results for Mi = 400,600 and 800 GeV. The general trend seen in Fig. 11 
is that increasing |5rj.| implies a lower instability scale Auv. When the instability scale is very low, however - 
lower than 10 TeV in these examples - we note that some of the lines in Fig. 11 curve backwards and indicate 


a slightly larger Auv for larger value of dr,.. This is the threshold effect that was mentioned in the previous 
paragraph; again, this behavior is only apparent where Auv is very low in the first place. 

As already seen in the case of colored fermions, when the Yukawa couplings defined at the vector-like 
fermion threshold scale are sufficiently large, then their subsequent RGE leads to a break down of our pertur¬ 
bative calculation on a scale smaller than the nominal vacuum instability scale. When this happens, we repeat 
our procedure of the previous section and define Auv as the scale at which either | F| or | F'^ | reaches a magni¬ 
tude This is visible in the left panel of Fig.[ll|for Mi = 400 and 600 GeV and dr,. > 0.1 and 0.2 respectively. 




The dashed lines show the results when modifying the perturbativity benchmark 4n/V4D by 20%. 


4 Conclusions 

In this paper we studied constraints on theories beyond the Standard Model, containing new physics fermions 
but no new scalar or vector-boson particles up to a very high scale. The basic constraint we discuss is due to the 
vacuum stability of the Higgs effective potential. In the presence of new fermions coupled to the SM through 
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Figure 11: Left panel: the maximum Ryf coupling deviation vs. the corresponding instability scale Auv, obtained by 
adding vector-like leptons with charged states carrying Q = 1. Purple, red, and orange lines correspond to fixing Mi = 
200,400 and 600 Gey respectively. The dashed lines correspond to modifying the perturbativity constraint by 20%. Right 
panel: same as on the left, but for Q — 2 and Mi = 400,600 and 800 GeV 


Yukawa interactions, the renormalization group evolution of the Higgs quartic coupling is negative definite, 
leading to vacuum instability. 

We surveyed Higgs couplings to the bottom and top quarks, to the tau lepton, and to gluons and photons, 
focusing on one particular coupling at a time. Our results imply that measuring a deviation at the level of 10% 
or so in any of these couplings (combined with null detection of Z coupling deviations in the Hbb or Htt case) 
would suggest a low cut-off scale Auv for pure fermion models, below about 100 TeV, where new bosonic states 
must stabilize the Higgs potential. While we did not investigate here the question of the actual mass scale for 
these bosons, the decoupling theorem leads us to expect that these states should occur not far above Auv- 
For large deviations [still at the ^(10%) level], the vacuum instability occurs essentially immediately above the 
fermion mass threshold, implying that the pure fermion model is not self consistent even as an effective theory. 

Our main numerical results are summarized below. To derive these results, that are expressed as consis¬ 
tency relations between Higgs coupling deviations and the vacuum stability cut-off Auv for pure fermion mod¬ 
els, we have judiciously tuned the parameters of the models in order to obtain conservative estimates, namely, 
maximal Auv per given Higgs coupling effect. For more generic parameter configurations we would expect the 
vacuum instability scale to occur earlier than these conservative estimates. 

• Hbb and Htt 

- For pure fermion models, vacuum stability constraints imply that measuring an Hbb deviation 
\5rh\ ^ 0.2 would require a corresponding deviation in theZbb couplings at the permille level, well 
within the expected sensitivity of future precision experiments with ILC-like capabilities. Finding a 
large Hbb deviation without an accompanying deviation in Zbb can therefore be used to rule out 
pure fermion models. The result for the Htt coupling is similar. 

- Currently available Z-pole data from LEP imply that partial vector-like fermion representations - 
models in which new fermions mix with either the SM left-handed quark, or the SM right-handed 
quark, but not both - cannot induce \5r}j\ > 0.1. This result is independent of vacuum stability 
arguments. The result for the Htt coupling is similar. 
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• Htt 


- Vacuum stability constraints imply an upper bound on tbe possible Htt coupling deviation in pure 
fermion models that we estimate, conservatively, as 15 rj | < 0.25 for Auv >100 TeV 

- In contrast to tbe Hbb and Htt examples, current precision electroweak data put only mild con¬ 
straints on tbe effective Ztt coupling, meaning that a partial vector-like fermion representation, 
mixing with either left- or rigbt-banded SM quarks but not both, could still induce a sizable Htt 
deviation. In sucb case, however, there should be a corresponding comparable deviation in Ztt. 


• HGG 

- Vacuum stability constraints on HGG deviations induced by pure fermion models depend crucially 
on tbe experimental mass limits on tbe new states. Allowing for color triplet fermions as light as 
500 GeV we find that imposing that the vacuum instability scale Auv > 100 TeV implies Ifird < 0.2. 
If the new states are more massive than 1 Tey imposing Auv > 100 TeV implies Idrc | < 0.05. 

- High color representations are also constrained. Allowing color octets or sextets as light as 1.2 TeV, 
we find that Auv > 100 TeV implies 5 re < 0.2, with somewhat weaker constraints on the magnitude 
of negative d rg. If the new states are more massive than 2 Tey the same Auv bound implies | d re | < 
0 . 1 . 


• Hyy 

- Considering new fermions with electromagnetic charge |Q| < 1, and allowing for states as light as 
200 Gey we find that Auv > 100 TeV implies —0.3 < dr^ < 0.15. If the new states are more massive 
than 400 GeV, the same Auv bound implies -0.1 < dr^ < 0.05. 

- Considering new fermions with exotic charge |Q| = 2, and allowing for states as light as 400 GeV, 
we find that Auv > 100 TeV implies 5rj-< 0.2, with somewhat weaker constraints on the magnitude 
of negative dr^-. If the new states are more massive than 800 Gey the same Auv bound implies 
Idr^l <0.1. 

Finally, a few comments are in order. The experimentally accessible signal strength (production cross sec¬ 
tion times decay branching fraction) in various analysis channels depends on more than one underlying ef¬ 
fective Higgs coupling. For example, the HGG coupling affects the signal strength in H —> yy measurements 
by modifying the gluon fusion production cross section, and the Hbb coupling affects all other signal strength 
measurements through its effect on the total width and thus on the respective branching fractions. Glearly, 
moreover, a realistic new physics scenario may involve true deviations in a number of different Higgs cou¬ 
plings, adding to the complexity. In this paper we chose to ignore this complication in the interpretation of 
Higgs data, assuming simply that the various degeneracies can be resolved sufficiently well. Our theoretical 
vacuum stability constraints apply therefore to each underlying coupling individually, and should be cast into 
constraints on signal strengths once new data becomes accessible, during Run-II of the LHG or with future 
colliders. 

Our calculation of the vacuum stability constraints relied on perturbation theory. In some models and 
some corners of the parameter space, the renormalization of the Yukawa couplings leads to a breakdown of the 
perturbative calculation on scales below the naively deduced vacuum Instability scale. We indicated where this 
happens in the relevant sections of the HGG and Hyy analyses. In practice, when this happens, the relevant 
scales are very low, in the ballpark of 10 TeV. 

We did not investigate the possibility of adding several different (non-minimal) vector-like fermion repre¬ 
sentations, and using them, e.g., to add a multiplicity factor in the HGG or Hyy analyses or to cancel some 
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of the correlations between Higgs and Z observables in the Hbb or Htt analyses. Some considerations along 
these lines can be found in . Since adding more fermions and more Yukawa couplings tends to make the in¬ 
stability problem more severe, we expect that, without significant accidental cancelations, such non-minimal 
models would result in vacuum stability constraints that are comparable to or stronger than those we derived 
here for minimal models. 

Our account of the experimental constraints on the various vector-like fermion models was partial, omit¬ 
ting potentially important constraints from, e.g., electroweak oblique parameters (see e.g. [^|^|^|^) and 
precision flavor data, including these additional constraints can only make our results stronger. 

Lastly we comment on deviations to the HWW and HZZ couplings, in a pure fermion model, contribu¬ 
tions to the effective HWW and HZZ couplings arise at the loop level. In contrast to the HGG and Hyy cases, 
however, where we have analyzed the analogous loop contributions, here to estimate the Higgs coupling de¬ 
viation one should compare a tree level coupling with a loop-suppressed effect. This means that the vacuum 
stability constraints on pure fermion models are strong, sufficiently strong to imply that observing modifica¬ 
tions to the HZZ or HWW couplings with the precision expected at the LHC would most likely rule out any 
pure fermion model interpretation. 
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A The effective potential 

In this section we describe our calculation of the Higgs effective potential in the presence of the canonical 
vector-like fermion representation 

j(D,1)q, j''(D,1)-q, (A.l) 

adding to the SM the potential 

'Ynp = YH'^ipx+ Y‘^H^eip‘^ x'^+ + 

For simplicity, we do not include mixing of the heavy vector-like fermions with SM fermions. While this mixing 
is important for Higgs coupling deviations, the vacuum stability constraints that we are concerned with in 
this section are easier to derive without it, and provide constraints that are easily generalized to the mixing 
case in the phenomenologically interesting regimes. We also restrict our analysis to the case where all of the 
parameters m^, m^, Y, and Y‘^ are real. 

We use the MS one-loop Coleman-Weinberg (CW) potential in the Landau gauge improved by two loop 
RGE (see e.g. fl4}[T7] ). Our problem involves two distinct physical scales - the SM electroweak scale, that we fix 
here at Mf ~ 173 GeV; and the vector-like fermion mass scale, denoted by M, that we define as 

M=.^(m2+mp/2. (A.3) 

In most of our analysis, M » Mf. We thus need to implement explicitly the decoupling of heavy states in the 
calculation of the effective potential, as we now explain. 
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The one-loop effective potential is a function of the classical Higgs field he and the MS renormalization 
scale pt. We implement the decoupling of heavy states hy defining two versions of the theory, one valid at 
jx <M and the other valid at pt > M, and matching between these theories at pt = M. In both of the regimes, 
pt < M and pt > M, we define 


V=Vo + Vi, 


(A.4) 


where the tree level piece is given by 


Vo=n[t)- 




For the one-loop contribution, on scales po M we use 






(A.5) 


(A. 6 ) 


where the sum on i includes W,Z,t, the Higgs h and Goldstone bosons and the heavy vector-like fermions. 

The MS constants Ci are equal to 5/6 for vector bosons and 3/2 for fermions and scalars. The effective num¬ 
ber of degrees of freedom m are nz = 3, nw = 6 , nt = —12 and n = —4D for each new 

fermion mass eigenstate transforming under SU(3)c in representation D. All couplings are running as function 
of the RGB parameter t = log [fi/Mt] with two loop RGB that we derive using the package provided by |[^ . The 

renormalized Higgs field is h[t) = ^[t)hc with ^[t) = \ where y is the Higgs anomalous dimension. 

On scales pt < M, we replace Bq. fA. 6 | by 




teSM 


6471^ 


mf(h,f)log 


m?(h, f) 


cqH . 


(A.7) 


Here, the sum on i runs over the SM states, excluding the heavy vector-like fermions. The couplings run with 
beta functions and anomalous dimensions obtained with the SM RGB equations. The Ceh® term contains the 
leading non-renormalizable dimension six operator induced by integrating out the heavy states. To obtain it, 
we expand the contribution of the heavy fermions to Bq. fA. 6 | in powers of the Higgs field h, fixing jx = M, 
and reading off the coefficient of h®. The full expression for ce is not particularly illuminating: in the case 
= rriip = M and = F, for example, it is given by Ce = DY^/{960n^M^). We have verified that adding 
higher order operators (i.e. csh®) does not affect our results. _ 


While we use the same symbols for the parameters A, and H in Bq. (A.5 , above and below the matching 
scale pt = M, their numerical values are different due to the one-loop threshold corrections. Writing A(pt = 
M-I-) = A(pt = M-) -I- AA, m^(pt = M+) = m^(pt = M-) -I- Am^, and f2(pf = M+) = f2(pt = M-) -I- Af2, we have 
AA = —4 c4 , Am^ = 2 c 2 , and AH = —co, where the corrections c„ are given again by the coefficient of h" on 
expanding the vector-like fermion contribution to Bq. (A .6 in powers of h. In the case = niip = M and 
F*^ = F, for example, we have Cq = 3DM‘^/[47iy, Cz = 2DY^M^/{4n)^, and C 4 = —DY^/{l2n^). 

We fix the RGB initial conditions for the SM couplings at the scale pt = Mf. While our calculation only 
requires one-loop matching, for convenience we use the two-loop values for the SM parameters at pt = Mf as 
given in Ref. 169 . For and A(t) we include the non-renormalizable Ce contribution through the correction 


A( 0 ) 

m 2 ( 0 ) 


A(SM|(o,+ ^ 


■f fdAVA 

I Sh 

fdAVA 1 

I 3h 2 


. Jh=v. 

52Ali^ 

Jh=v’ 


(A. 8 ) 

(A.9) 
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with AVi = ceh^ and where and are taken from Ref. . Lastly, for the vacuum energy f2, we 

set f2(0) = 0. This guarantees n(f) 0 up to logarithmic corrections at all t. 


The input values for the vector-like fermion parameters appearing in Eq. A.2 are taken to he defined at the 
scale /i = M. We note that, in principle, we should run the non-renormalizahle Ce contribution at /i < M using 
the SM anomalous dimension, as well as include insertions of Ce in the RGE for the SM couplings (5^ 60 . In 
practice, our interest in the effective potential mainly concerns scales /i > M, where vacuum instahility occurs, 
and where the non-renormalizahle operators are replaced hy the full one-loop contribution that we evolve 
properly using the two-loop RGE. The contribution of Ce is thus mainly in setting the initial conditions for the 
running of A. Here, however, the running of Ce would be a two-loop correction to the one-loop matching that 
our approximation requires, and so we omit it. Concerning the insertion of Ce in the running of A, using the 
results of Ref. [M ■ 60 we verified that including this effects leads to negligible corrections to our results. 

Now that we have the effective potential V[hc,n) defined for all /r, the final step in the calculation is the 
RGE improvement, that amounts to letting pt be a function of he in order to control large logs far in field space. 
Wesetpt= he- 

Einally we note that in some of the analyses in the body of the paper, in parameter regions allowing for 
very light fermions, the vector-like mass scale M comes out lower than the SM top quark mass Mf (in practice, 
this only happened for the Mi = 200 GeV curve in the left panel of Eig. 


11 


and only for < 0). In this case, 
instead of the matching procedure described above, we used the full theory effective potential and RGE starting 
from Mt and setting the RGE initial conditions based on Eqs. i A .8 A.9 1 with A Vi replaced by the full log 
contribution of the new fermions. 


B EFT analysis of deviations in Z and Higgs couplings to SM fermions 

We examine the vector-like fermion representations that can modify the effective Higgs Yukawa couplings at 
tree level (for related analyses see, e.g. |^70||71| ). Before we introduce the new fermion fields, we first describe 
our conventions for the SM effective theory below the vector-like fermion scale. 

Our notation for the SM Higgs and fermion representations is iT(l, 2 ) i, q{3,2)i, d^(3, l)i , u^(3, 1)_2, /(l,2)i, 
e^(l, l)i. The first and second numbers in parenthesis denote the SU{3)c and SU{2)w representation, respec¬ 
tively, and the subscript denotes the hypercharge. The fermion sector of the SM Lagrangian is 

■^SM = iqo'nSs^q + 

-[ydH^qd‘^ + yuH^equ‘^ +yeHUe‘^ + cc^ (B.l) 

with €\2 = —621 = 1. In unitary gauge H = {0 h)^/^/2, with (fi) = i; ~ 246.22 GeV. 

Integrating out heavy fermions leads to effective Hff and Zff couplings that are modified compared to 
their SM values. The main effect is captured by considering non-renormalizahle operators (7^[7^ , 

A^eff = Ciffi + CC. (B.2) 

Following Ref. and adding operators that affect Yukawa couplings, we list the operators of interest as fol- 


21 












low^ 

6hi = iti'^^i^HlcT^l, = i , Ghe = iH^^^He‘^CT^e‘^ 

Ghq = iH^^i^Hgai^q, 6'^^ = i ■ {qa^^aq) 

Ghd = d‘^aij,d‘^, Ghu = 

GYd = \H\^H^qd‘^, GYu = \H\^HTequ‘=, GYe = \H\^HHe<^. (B.3) 

For fermion f, we define the modification to the Higgs coupling as 5rf = (^ghff — sff/j /sf/f We find 

^Note that our definitions for the operators &i,d <Ghu< and 6ue differs hy a sign compared to the corresponding operators in Ref. . 
In addition, we define v = 246 GeV while Ref. works with f = 174 GeV 

Table 1: Vector-like fermion representations that can mix with SM quarks 



representation 

Ynp 

DI 

Q(3,2)i, Q"(3,2)_i 
D(3,1 )!i, D"(3,1)I 

3 3 

YQdcWQd'^ + YqocH^qD^^ + YqdcWQD^^ + YqcoHTeQ^^D 
+MqQ^€Q‘^ + MdDD‘^ + cc 

DII 

Q'[3,2)_5, Q'"(3,2)| 
D(3,l)_i, D"(3,l)i 

3 3 

YQ,dcHT€Q'd‘^+ Yqo^WqD^^ + Yq, eQ> D<= + YqUoH^Q'‘^D 

-FMg/Q'^eQ'^ -FMflDD'^ + cc 

Dill 

Q(3,2)i, Q"(3,2)_i 
D'[3,3) 1, D'^[3,3)i 

3 3 

YQdcWQd^ + Y^qj^cWaq-D'‘=+YQj^^,H^(7Q-D'^ + YQciyHT€aQ‘=-D' 
-FMgQ^eQ^ + MjyD' ■ + cc 

DIV 

Q'[3,2)_|, Q'‘^(3,2)| 
D'(3,3) 1, D'^(3,3)i 

3 3 

YQ>dcH'^€Q'd‘^ + Y^jj'.H^aq ■ D’^ + Y^j^'M'^eaQ' ■ D’<^ + YQ>,jyH^aQ'‘^ ■ D' 
+Mq'Q^'^€(^‘^ + Md'D' -D'^^ + cc 

UI 

Q(3,2)i, Q"(3,2)_i 
1/(3,1)L m(3,l)_2 

3 3 

YqucH^ eQu^ +Yqu-HT eqW +YqucH^ eQU^^ + YqcuH^Q^U 
+MqQ'^€Q‘= -F MuUW -F cc 

UII 

Q"(3,2)7, Q"‘^(3,2)_7 
U{3,l)l U‘={3,1)_2 

3 3 

FQ«^.//tQ"u‘^ + F^yc/Z^et/m-F FQ«y,//tQ"m + Fg/ZcyH^eQ""!/ 
-FMQ«Q"^eQ"^ -F MuUW -Fee 

mil 

Q(3,2)i, Q"(3,2)_i 
U'{3,3)2, U'‘^{%3)\ 

3 3 

YqucH^€Qu‘^ + Y^jy.HTeaq ■ t/" + Y^^M^ecrQ ■ + Yqcu'H^(tQ‘^ • U' 

+MQQTeQ‘^ + Mu'U'-U'‘^ + cc 

UIV 

Q"(3,2)7, Q"‘^(3,2)_7 
U'{3,3)2, I/‘'(3,3) 2 

3 3 

Yq« + Y^^y,HTecTq ■U'‘= + YQ>^^,WaQ" ■ U'<^+ Y^^,^j,H^€a(^'<^-U' 

-FMQ«Q"^eQ"'^ -F Mu' U'■ if‘^ + cc 

DV 

Q'(3,2)_|, Q'‘^(3,2)| 
D"(3,3) 4, D"^(3,3)4 

3 3 

Yq'dcH’^eQ'd'^ + Y^j^»cH^aQ' ■ D"‘^ + Yq'cj^"H'^ eaQ''^ ■ D" 
-FMg/Q'^eQ'^ -F Md"D'' ■ D"‘= + cc 

uv 

Q"(3,2)7, Q"‘^(3,2)_7 

U"{3,3)1, l/'^(3,3) 5 
_2_2_ 

YQ",^cH^(^'u‘^ + YQ"^'cH^ea<:/' ■ U"<^+ Yq"cu "• U" 

-FMQ//Q"^eQ"^ -F Mu" U" -if'^^ + cc 
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(see also fsi] ) 




y ^hd^ f 5gAd — 2 * 

~ ^^hq ^iiq ^hu^ > ^gAu — ^ (^^hq ^j^q ^hu^ > 

\ ,.? r . \ 


?Vd— 2 i^^hq + c'^^ 

-£^Cy., 5gy.=-^^ 

5^6 — ^SVe — ^ ^/z/ ^he^ } ^gAe — ^ {^^hl ^ ^he^ } 

^Svv — Y ^/i/) ' ^Sav — Y 

In addition, for any /, the ( 


,_, r r Z V 

5 couplingy/ of Eq. ^B.l| is constrained by 

w~f'-T 


^hl 


.)■ 


(B.4) 


(B.5) 


Table 2: Contributions to the non-renormalizable operators involving SM quarks, listed in Eq. B.3 



Effective operators 

DI 

_ „ _ \YqD‘:\^ _yd\YqD‘:\^ . JdligdCp YqC q)YQac 

Chd— 2 \MqP’ ^hq~^hq— 4|MdF ’ — 2\MdP 2\Mq\^ MqMd 

DII 

_ liQd'^P . _ \YqD<^\^ _yd\YqD‘:? ydliQdCp . YqC j^Yq^c Y^o'^ 

Ciid— 2 |MqP’ 4|MdF ’ — 2|MdP 2|MqP ' MqMd 

Dill 

\YQdc\^ ^ , 3|V,ocp yd\YqDc\^ ydliQdCp YQCqjYQacY^oC 

2|Mq|2’ ‘'^'7 4|Md|2 ’ 2|Md|2 2|MqP MqMd 

DIV 

^ _|iQd'^P ov _ ^ _3 |V,d<:F _ydlV'l" 1 ydliQd'l" 1 YqcdYq^cY^dc 

2|MqP’ 4|MdP ’ 2|MdP 2|Mq|2 ' MqMd 

UI 

_\Yqu‘:P _yu|i'^[;cp yaliguCp YqCuYu^cY^uc 

C/iU— 2 |Mq| 2’ C/jq— 4 |Mj,| 2> Ciu— 2|Mt/|2 ' 2 |Mq|2 ' MqMu 

UII 

_ IYqucP _ yu|5ij[/c|2 ^ yiiliijuCp YQCuYuqcY^uc 

Chu— 2|Mq|2’ C;2(?— 4|M[,|2’ 2|M[,|2 ' 2|MqP MqMu 

UIII 

_ I^QucF o / _ _3|V,[;c|2 _ yaligucp yuli'^t/cp Yqc ijYu^c Y^uc 

C/JU — 2|Mq| 2’ ^ “ 4|Mt,|2 ’ 2 |Mq| 2 2|Mt/p MqMu 

UIV 

_ I^Qu'^P Of,/ _ f, _ 3|i;j[7‘;F _yu|iQiiCp YulYquCp , YqC uYu^c Y^ui: 

Chu— 2|Mq|2’ ~ ^^<7 ~ 4|M[7|2 ’ 2|Mq|2 2|M[,|2 ' MqMu 


23 



























































Eqs. I B.4 B.5 are valid to [V^ v^/M^), where M represents the heavy vector-like fermion mass scale and Y is 
a Yukawa coupling. 

Vector-like fermion representations that can mix with SM quarks are given in Tab. Integrating out the 
heavy fields we collect contributions to effective operators of interest in Tab. This is further summarized in 


Table 3: EFT contributions to Higgs and Z couplings to SM quarks. The complex phase (j) has a similar structure in all 
cases; for example, in model DI it is given by ^ = arg . 
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Table 4: Vector-like fermion representations that can mix with the SM leptons 



representation 

Ynp 

LI 

L[l,2)_i, L"(l,2)i 

YLecWLe‘^ -[-YiecH^IE<= + Yle^WLE‘^ + YLCEH'^eHE 
+MlL'^€L‘^ + MeEE‘^ + CC 

LII 

L'(l,2)_|, L'^(l,2)| 
E(l,l)i, 

YL>e'=H'^eL'e‘^ + YiEcHUE‘^ + YL>EcHT€L'E‘= + Y^>cEH'^L'‘=E 
YMEL'TeL'^ + MeEE‘= + cc 

LIII 

L[l,2)_i, L"(l,2)i 

F(l,3)i, E'"(l,3)i 

YLecH^Le‘^ + -E Y^j/cH^ctL-E'‘^ + Yece'H'^ eaH ■ E' 

YMeLT' eU + Me'E' ■E'<= + cc 

LIV 

Z,'(1,2)_3, L'‘^(1,2)3 

F(l,3)-i, E'"(l,3h 

YE^cH^eL'e‘= + Y^j,uH^<Jl ■ E'‘= + Y^,j,>cH'^e(jL' ■ E'‘= + • E' 

+Me L''^eL'‘= + ME'H-E'‘= + cc 

LV 

L[l,2)_i, L‘^(l,2)i 
E"(l,3)o, E"^(l,3)o 

YiecWLe'^ + Y^j/cH'^eal ■ E"<= + Y^ce"H^(tL‘^ ■ E" -E eaL-E"‘^ 

+MlL'^€U + Me«E" ■ E"‘= + cc 


Table 5: Contributions to the non-renormalizable operators involving SM leptons, listed in Eq. B.3 



Effective operators 

LI 

,,, _ -r, Ihfid" ,, _ Y„c\Yiec\^ Yi,c\Yi,,c\^ Yi,ceYl,cY,ec 

2|MiE’ ^hl 4|Mi;p’ 2\Me\^ 2|MlE MlMe 

LII 

_ Ihfid" „ _ h.dhfid" , Yi,c\Ye,c\^ , YeceYe.cYiec 

2|MiE’ ^hl 4\Me\^’ 2\MeP 2|MlP ' MlMe 

LIII 

_ _ „ _Y„c\Y,Ecp YeceYe^cYiec , 

2|MiE’ 4\Me\^ ’ 2\Me\^ MeMe 2|Mip 

LIV 

_ Ifi.d" o„/ _ „ _Yi,c\Y,ECp YeCeYe^cYieC Yi,c\Ye,c\^ 

2|MiE’ A\Me\^ ’ 2\Me\^ MeMe 2\Me\^ 

LV 

_ \YEeCp _ 3|y;B.|2 _ Yi,c\Ye,cP Y„c\YiECp n YeC eYe,c Y, eC 

2|MlP’ ''hi 4\Me\^ ’ 2|MlP ' |MbP MeMe 


terms of relations between the Hff and Zff effective couplings in Tab.|^ Note that the expressions in Tab.|^ 
include an expansion in Idg^l and |5gv|, and are valid to leading order in the 5g’s. 

Vector-like fermion representations that can mix with SM leptons are given in Tab. Integrating out the 
heavy fields, we further summarize the contributions of the effective operators in Tab.|^and the modifications 
to the Hff and Zff effective couplings obtained via the EFT analysis in Tab.[^ again valid to leading order in 
|5gyi,vl- These results are used in Sec.|^ where we also comment on cases in which some of the expressions 
derived here using EFT remain valid even for 
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Table 6: EFT contributions to Higgs and Z couplings to SM leptons. The complex phase 0 has a similar structure in all 
cases; for example, in model LI it is given by 0 = arg {^YiceYi e‘ Yie‘y*M*j^M*^'j. 



Coupling deviations 

LI 

Ore~ ZOgAe+ '\]\^Sve ^ ° SAe 2 \2\Me?- 2\Ml\^ ) ’ ^ SVe 2 l^2|Mi;|2 2\Ml\'^ ) 

LII 

Ore~ 2ogVe _f2me ^Sab]^ ’ ^ SAb 2 1^2|MbP 2|MlPJ’ ^ SVb 2 2|MiJp 2|MlP J 

LIII 

Ore~ 2ogVe y2me yr^VB ^ Sab\^ ’ ^ 8Ab 2 2|Mi7|2 2 |MlP J ’ ^ 8Ve 2 \2\Me\^ 2\Me\^ j 

5gAv = 5gVv='’2 

LIV 

Ore~ 20gAB+ tpirm yr^VB ’ hg/le 2 i^2|MBp 2|MlPJ’ ° 8Vb 2 y2\ME?- 2|MlPJ 

LV 

Xr '-''’XiT 1 /|XrT2 X rr2 D'h X rr _ T Ihfirp . IlLeCpA ^ _ v'^ ( Ihcrp IrteCpA 

Orer^2ogVB+ y2me yr^VB ^8 ab\^ ’ OgAB 2 IMbP "*■ 2|Ml|2 J ’ ^8Vb 2 \MeP 2|MlPJ 

5gAv = 5gvv= 


We comment that the derivative operators in Eq. jB.3| , for which we have highlighted the effect on Zff 
couplings, also enter into Higgs ti associated production and h bb,TT decays. However, in contrast to the 
modified Yukawa couplings [parametrized by the 5 rf terms in Eq. l |B.4| ] , the contribution due to the derivative 
operators does not interfere with the leading SM contribution to the Higgs decay or production matrix element. 
Thus, for the Htt, Hbb, and Htt couplings of interest to us in this paper, their contribution to the respective 
signal strength is suppressed in comparison to the 5 rf terms, and we neglect it here and in the main text. 


C Collider constraints 


We discuss here the collider bounds on vector-like quarks and leptons. We also include limits on states with 
exotic color or electromagnetic charge assignments that can be relevant for the Nyy and HGG analyses. 

Unless specified otherwise, we list constraints on a single Dirac fermion at 95% C.L., taking into account 
the color multiplicity, but ignoring additional constraints due to e.g. isospin multiplicity. Eor the vector-like 
fermion models at hand, these constraints are conservative since they ignore the contributions due to some of 
the states. Thus our bounds are weaker, for example, than those derived in for vector-like leptons. 


We discuss the constraints on states that are either stable or decay promptly on collider time scales, sepa¬ 
rating the discussion for colored and color singlet fermions. We do not discuss the constraints on states that 
decay non-promptly but yet within the tracker volume. We comment however that for ct > 0.1 cm the con- 
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straints on such non-prompt decays are typically comparable to or stronger than the constraints obtained in 
either the stable or prompt cases (see, e.g. (77|[7^ ). 

We stress that the broad summary of constraints presented here, that applies to numerous different vector¬ 
like fermion models, is only intended to provide a rough estimate for the mass scales of such fermions still 
allowed after Run-I of the LHC. For the purpose of the current paper we find these rough estimates sufficient; 
a more careful analysis will be motivated in the case that Higgs coupling deviations end up being discovered at 
the LHC. 


C.l Colored particles 


As discussed in Section |3.1| we are interested only in states transforming as 3, 6 or 8 ofSf/(3)c. For the vector-like 
models we deal with here, higher color representations induce a Landau pole in Ug on scales close to the vector¬ 
like fermion masses. To obtain the mass bounds quoted in the following, for the color 3 we use the t' cross 
section computed in [7^ using HATHQR and for the 8 we compute it at NLO-FNLL using NLLfast |81 
The mass bounds for the 6 and 8 representations are similar 


C.1.1 Stable states 

If we are interested only in iTGG or iTyy deviations we can imagine that the new fermions possess a conserved 
quantum number that prevents them from decaying into SM states. In this case we are looking for stable 
colored particles that travel inside the detector and can form charged bound states. The strongest bounds 
on heavy stable charged particles are set by the CMS search that exploits mainly their long time of flight and 
their anomalous energy loss per unit length [^. The behavior of heavy colored particles inside the detector is 
subject to large uncertainties, so we quote a range for the mass exclusion given by the different hadronization 


models considered by the collaboration ■ 90 . The bounds we quote here refer to a combination of the CMS 
tracker and time of flight (TOF) analyses. Similar bounds are obtained from the tracker analysis alone. 

The CMS collaboration considered colored particles with Q = 0 and Q = 2/3. If we use the cross section 
limit obtained by the collaboration on gluino pair production, we find the mass bounds m 3 > 1050 ±30 GeV 
and mg > 1295 ± 35 GeV for a single Dirac fermion and m 3 > 1165 ± 35 GeV and mg > 1390 ± 35 GeV for three 
Dirac fermions (corresponding to a vector-like doublet and a singlet with the same mass). The bounds are the 
same for Q = 2/3 and Q = 0 states. We expect the limit on the 6 to be similar to that on the 8. This procedure is 
partially justified by the fact that the cross section limit for stops (i.e. Nc = 3 and Q = 2/3) converges to the one 
for gluinos (Q = 0 and Nc = 8 ) at high masses. However it can overestimate the bound for particles with Q = 0 
and Nc = 3,6 and a more thorough collider study is needed to go beyond this very rough estimate. We expect 
similar bounds for other electromagnetic charges, including e.g. Q = 1/3. 


C. 1.2 States decaying promptly to first and second generation quarks 

Wq, Zq and hq final states Vector-like quarks with mass mixing with the first or second generation SM quarks 
can decay to W,Z or h plus one jet. This situation could be relevant for the HGG andiTyy analyses, where mix¬ 
ing with third generation quarks is not guaranteed. Bounds on pair produced heavy quarks for these final states 
were examined in The corresponding experimental analyses where not updated since the 7 TeV run, and 
the leading explicit constraint we find is m^' > 350 GeV from ATLAS |^, where the search was performed on 
WW±2 jets final states and the bound assumes BR(^' —> Wq) = 1. For three Dirac copies of q', the ATLAS 
search would exclude m^/ > 420 GeV ATLAS also reported a multijet search for RPV gluinos |j^ that is a count¬ 
ing experiment, does not exploit the shape of the jet invariant masses and can be reinterpreted in the context 
of our models. A possible caveat is that in most cases, ATLAS finds that requiring a seventh jet improves the 
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bounds. Our signal, that includes WWjj, ZZj j and hhjj final states, would has a different probability of ra¬ 
diating an extra jet with respect to a pair of gluinos, but we ignore this correction for the purpose of the current 
rough estimate. In the case BR(^' ^ Wq) = 1, no fi-quarks are produced and the RPV gluino search yields a 
weaker bound than the one we obtain from 1^ discussed above. Instead for BR(( 7 ' —> hq) = 1 we get for three 


Dirac fermions niq' > 700 GeV while for BR(q' —> Zq) = 1, niq' > 450 GeV. The ZZjj and hhjj final states 
could also be constrained by the 8 TeV CMS multilepton search performed with an integrated luminosity of 
19.5 fb“^ 194 . For BR(q'' ^ hq) = 1 we find that the mass bounds for three Dirac fermions is m^/ > 540 GeV, 
which is weaker than the multijet bound quoted above, while for BR(q' —> Zq) = 1, mq> > 650 GeV; stronger 
than the multijet bound. 


Dijet decays Both the 3 and the 6 could couple to an SM quark and a gluon via nonrenormalizahle magnetic 
interactions, allowing dijet decays. However we find somewhat unlikely for the dijet mode to be the dominant 
one, as producing the nonrenormalizahle magnetic interaction requires flavor changing scalar or boson loop 
diagrams. If the operator is formed with a W boson loop (relevant for the 3), then we expect the irreducible 
q' —> Wq to dominate and lead to the hound discussed in the previous paragraph. For the 6 , new bosons with 
masses larger than Auv must be involved and we expect a tree level three-jet decay (discussed below) with 
branching fraction comparable or larger than the dijet mode. However from the the CMS search for pairs of 
dijet resonances we can make a rough estimate m 3 > 500 GeV, valid for a single Dirac fermion, with a 
somewhat stronger limit on the 6 . 


Three jet, and higher jet multiplicity decays Limits on this final state were set by the ATLAS multijet search 
discussed above |^, with somewhat weaker bounds reported by CMS 1^. For a single Dirac 8 we find ms > 
1.1 TeV; while for three copies of equal mass we obtain a limit stronger than ms > 1.2 TeV. The latter bound is 
conservative, since what prevented us from quoting a higher number is simply the fact that the ATLAS exclusion 
plot extends only up to 1.2 TeV The constraints on the 6 representation are the same. For the 3 representation, 
we find essentially no bound from . However, our new fermions necessarily include electromagnetically 
charged states, leading to irreducible three jet plus weak gauge boson (or Higgs) decays which translate with 
ff{l) branching fraction into five jet modes. From the 10 jet bin of the general purpose analysis in j^, we find 
m 3 > 850 GeV or so. In addition to the five jet channel, leptonic modes [with an ^(10%) branching fraction] 
suggest even stronger bounds of order a TeV 


C. 1.3 States decaying promptly to third generation quarks 


The discussion here is aimed to address vector-like quarks mixing with the third generation SM quarks, relevant 
to the Hbb and Htt coupling analyses. The constraints in this case are generally stronger than those discussed 
above for decays to W,Z,h + q, involving only first or second generation quarks. We focus on the CMS searches 
for pair production. CMS results are presented as a function of branching ratios and can be directly read off 
the tables in ^ 9^ . ATLAS bounds are comparable [100 jl03| . 

In the case of a t', the weakest bound is rrit' > 700 GeV] obtained for BR(t' —> Wb) 1 (^. Adding a finite 
BR(t' —»Zt) or BR(t'—» ht) makes the hound a few tens of GeV stronger. For a b' the situation is reversed, since 
top rich final states are easier to distinguish from the background. The most constraining case corresponds 
to BR(h' —> Wf) 1 and mh> > 730 GeV] while the weakest bound reads mf > 600 GeV, obtained for BR(h' —> 
Wt) = 0 andBR(h'—> hb) f^BR{b' Zb) |[^[^ . 

Even if the lightest state has an exotic charge, the bounds do not change dramatically. For a quark with 
charge 5/3, we have BR(t®/^ —^ Wt) = 1 and m^s/s > 800 GeV |l04 
the limit is the same as in the t' case discussed above. 


A charge 4/3 quark will decay to Wb, and 
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C.2 Leptons 


We focus on pair production pp —» LiLi where Li stands for the lightest new lepton, that in our minimal 
models is always electromagnetically charged if Higgs coupling deviations are to be induced. Our limits on this 
single state are less stringent than those set on complete models by (M and should be thought of as 


conservative estimates. We consider separately the cases in which Li is mostly an SU{2)w singlet, doublet or 
triplet. 

We compute the relevant cross sections at QCD NLO with Pospino |105| decoupling all supersymmetric 
partners except for Winos (Higgsinos) when considering SU(2) w triplets (doublets). The singlet Drell-Yan pro¬ 
duction cross section was computed at LO in using PYTHIA v6.426 1 106] and CTEQ6L1 PDFs |107 . For 
doubly charged leptons we generate the models in UFO format |108 using FeynRules 1109 110) and compute 
the cross section at LO with MadgraphS |lll) . 

We first discuss the case in which the lightest charged lepton is stable, relevant to the Hjy analysis, and 
then consider the case in which it decays to the SM states. For the latter possibility, the decay channels of Li 
are Wv,Zt(£), fiT[£) with branching fractions depending on its electroweak quantum numbers. For instance, 
if Li is mostly an SU(2)tv singlet or triplet and » m/j, it will decay to Wv,Zt(£), hT{^) with ratio 2:1:1. 
If Li is mostly an SU(2)w doublet, it will decay to Zt(£), fiT(£) with equal probability, but not to Wv, again for 
» m}i. In the following we use these ratios of decay probabilities to set our bounds. Note however that 
near kinematical thresholds there can be important differences. In particular, if Wv is the only allowed decay 
mode, LHC constraints essentially vanish, as discussed in Section C.2.4| 

While quoting mass limits on Q = 1 fermions we distinguish between singlets, doublets and triplets, refer¬ 
ring to F = 1 singlets, F= 1/2 doublets and F = 0 triplets. Bounds on F = 3/2 doublets can be inferred from 
the numbers presented here by increasing the F = 1/2 doublets cross section by 50%. Q = 1 fermions from 
F = 1 triplets instead have the same production cross section as Q = 1 singlets. 

In general we expect multiple leptons in the final state and we find that the strongest constraints at the mo¬ 
ment come from the 8 TeV CMS multilepton search performed with an integrated luminosity of 19.5 fb“^ . 
We compute the asymptotic CLs defined in the appendix of in each individual signal region of the search 
to constrain cr x e where cr is the production cross section and e includes acceptance, trigger and identification 
efficiencies, efficiencies of the kinematical cuts and the branching ratio of the vector-like lepton. In evaluating 
e, we computed the branching ratio for each signal region and assumed a flat 70% efficiency times acceptance 
for electrons and muons, and an hadronic tau identification efficiency of 35%, based on 1112 ■ 115 . We also 


assume a flat 70% efficiency for identifying a fi-quark We always quote the bound from the most con¬ 
straining search region. Clearly, stronger bounds could be derived in a more refined analysis by combining the 
results in different regions. 


C.2.1 Stable states 

The same CMS search relevant to stable colored states also sets the strongest bound on stable leptons |[^ . For 
Q = 1, the limits read mii > 574 GeV for singlets, > 670 GeV for doublets and > 800 GeV for triplets. 
For Q = 2 we find mi++ > 705 GeV for a singlet or doublet, and mi++ > 790 GeV for a triplet. 

C.2.2 Prompt decays to charged first and second generation leptons 

Considering Q = 1 states, a singlet L\ that only mixes with the first two generation leptons must have > 140 
GeV from the 4 lepton search region with missing Ej < 50 GeV; Hj < 200 GeV, one pair of leptons from a Z 
decay, no hadronic t and no b jets. When L\ is mostly a doublet the bound is mii > 260 GeV, from the 4 
lepton search region with missing Ej < 50 GeV, Ht < 200 GeV, two pairs of opposite-sign same-flavor leptons 
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where at least one pair originates from Z decays, no hadronic t and at least one b jet. For a triplet we find 
mil > 230 GeV from the same search region. For the case of doublet Q = 2 states decaying to W£, Ref. 
obtains mi++ > 460 GeV We expect the bound on singlets and triplets to be comparable. 

C.2.3 Prompt decays to charged third generation leptons 

Gonsidering Q = 1 states, in the singlet case, if Li only mixes with the t, the multilepton search bound on m^j is 
below 100 GeV, weaker than the LEP bound. If Li is mostly a doublet we have rriLi > GeV from the 4 lepton 
search region with missing Ej < 50 GeV, Hj < 200 GeV, one pair of leptons from aZ decay, no hadronic t and 
no b jets. If the lightest charged lepton is dominantly an SU(2) w triplet the mass bound is > 185 GeV from 

the same search region. For the case of doublet Q = 2 states decaying to MV I, Ref. obtains mi++ > 320 GeV. 

We expect the bound on singlets and triplets to be comparable. 

C.2.4 Prompt decays to a W and missing energy 

It is possible to add to our minimal models a SM singlet N (a sterile neutrino or a bino) mixing with the doublets 
L, via Yukawa interactions. This can improve the agreement with electroweak precision tests , though 
the additional Yukawa couplings imply somewhat more stringent vacuum stability constraints than would be 
obtained from the minimal models. In this setting the lightest charged Li lepton can decay predominantly to 
WN and N can be stable. If the decay is prompt, LEP still sets the most stringent constraint for small mass 
splittings between the neutral and charged fermion, giving rrii^ > 100 GeV for m^j — rriN > GeV. The 
ATLAS search for charginos is also relevant |116] , excluding masses up to 180 GeV for triplets, provided that 
m]v < 20 GeV If N is stable it could be searched for in mono-jet and mono-boson final states but current LHG 
bounds are not strong enough to constrain electroweak production. 

C.2.5 Prompt cascade decays of Q = 2 fermions 

In case of small mixing with the SM and largish mass difference between the new leptons, L++ —> L\ W+ can be 
the dominant decay mode. Assuming BR(L++ —> W+) = 1 and that Li only decays to the third generation SM 
leptons, we used the GMS multilepton search to obtain mi++ > (140 — 250) GeV for a singlet and mi++ > 
(180—300) GeV for a triplet. The lower end ofthe range is achieved when BR(L| —» Wv) = Iwhilethe upper end 
ofthe range is achieved when BR(L| — >Zt)= 1 orBR(L| —> fiT) = 1. Ifii only decays to the first two generation 
SM leptons, we have mi++ > (140 — 420) GeV for a singlet and mi++ > (180 — 490) GeV for a triplet, where the 
lower end of the range is achieved when BR(L| —> Wv) = 1 while the upper end of the range is achieved when 
BR(L+^fi£)=l. 
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